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In this work, the mechanism of Li insertion in 77'-Cu6Sn5 to form Li2CuSn is 
discussed in detail, based on both theoretical calculations and experimental re- 
sults. The mechanism is investigated by means of first principles calculations, 
with the full potential linearized augmented plane wave method, in combination 
with in situ X-ray diffraction experiments. The 77'-Cu6Sn5 structure, as well as 
its lithiated products, have been optimized and the electronic charge density cal- 
culated in order to study the change in bond character on lithiation. The average 
insertion voltage of the ?7'-Cu6Sn5-Li2CuSn transformation has been calculated 
to be 0.378V in good agreement with the experimental value. 
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INTRODUCTION 



Light weight and compact, hthium-ion batteries [J| are ideal energy storage devices for the use in appHances hke 
laptops, mobile phones and electric vehicles. But there are concerns about the safety of the lithium-ion batteries in 
their present form because, the most common material used as anodes in these batteries is graphite and its lithiated 
potential is very close to that of the lithium. The study of intermetallic insertion electrodes, as possible replacements 
for graphite, has recently attracted a lot of attention. Examples of such structures are InSb 0,0 and CuaSb H- 
The desirable characteristics for an ideal anode material include the capability of reversible Li intercalation/insertion 
while maintaining the stability of the structure. One of the materials that possess these qualities and has been found 
to be a possible alternative for graphite in lithium-ion batteries is r/'-CueSns 0, 0,0,13 ■ 

In situ X-ray diffraction (XRD) studies have shown that, during the electrochemical reaction of Ty'-CugSns with 
lithium, a Li2CuSn-type structure is formed 0,0 13- 'T^^ ideal reaction for this phase transformation can be described 
by: 



CugSns + lOLi — > SLiaCuSn -I- Cu. (1) 

Further lithiation of the Li2CuSn structure results in further extrusion of Cu from the structure, to finally form 
Li4.4Sn (the maximum lithiated Sn-phase). This reaction is described as: 



xLi + LiaCuSn — > Li2+xCui_ySn + yCu (0 < x < 2.4; < y < 1). (2) 

These reactions were recently supported by the in situ ^^^Sn-M6ssbauer spectroscopy measurements [3- When the 
lithiation process is interrupted, to only involve reaction (1) (above 0.2V vs. Li/Li"*"), the material shows a reversible 
and stable lithium insertion behaviour with a capacity of ^ 200mAh/g. If the material is fully lithiated to form 
Li4.4Sn, the cycling stability is dramatically affected and the ability to insert lithium reversibly gradually declines. 
The reason for this is believed to originate from the large structural rearrangements that occurs in reaction (2), along 
with the extrusion of Cu. The experimental data provides information on the main processes that occurs on lithiation 
of the Ty'-CugSns structure. The detailed mechanisms of reactions (1) and (2) are, however, unclear. For example the 
positions of the initial Li in the ry'-CueSns structure and the starting point of Cu extrusion from the structure are 
not known. Although there have been suggestions regarding the mechanism of the phase change from T^'-CugSns to 
Li2CuSn 0,13, there exists no experimental measurement or theoretical calculation to back these suggestions. In the 
present work we have performed ab-initio calculations to investigate this mechanism. The Ty'-CugSns structure has 
been optimized and the calculations performed to determine the average insertion voltage for the Cu6Sn5-Li2CuSn 
transformation. 

The paper is arranged in the following manner. Section II gives the details of the calculations and experiments. 
The results of the structural optimization and electric field gradients are presented in section III A. Section III B deals 
with the discussions about the charge densities and the nature of the bonds in ?7'-Cu6Sn5 and Li2CuSn compounds. 
In section III C the results and some ideas about the mechanism of lithiation of ?7'-Cu6Sn5 are presented. Finally 
section IV comprises of the conclusions of our work. 



METHODOLOGY 



Theoretical 



The total energy calculations are performed using the full potential linear augmented plane wave (FP-LAPW) 
method using the WIEN97 code 0]. The high-lying Cu and Sn semi-core d— states are treated using local orbitals. 
The calculations are performed within the the local density approximation (LDA). The scalar relativistic equations 
are used in a self-consistent scheme. These calculations are performed using 115 and 204 k points in the irreducible 
Brillouin zone (IBZ) for Ty'-CugSns and Li2CuSn respectively. The structural optimization for jy'-CueSns is performed 
using 4 k points in the IBZ. 

Ty'-CugSns has a monoclinic lattice with the space group C2/c (no. 15). The calculations are performed at the 
experimental lattice parameters [i3 (a=11.022A, b= 7.282A, c=9.827A). However, we have used the theoretically 
optimized atom positions given in table I for all the calculations. The conventional unit cell for ry'-CueSns is presented 
in figure 1. Structurally Li2CuSn exists in a cubic phase with the space group F43m (no. 216) The calculations 
are performed at the experimental lattice parameter of a=6.282A and the atomic positions are given in table I. 
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Experimental 

Ty'-CugSns was synthesised by a high-energy ball-miUing procedure, as described elsewhere 0. Electrodes were 
made by mixing 90 wt% CugSns powder with 5 wt% carbon black and 5 wt% EPDM rubber binder. The slurry was 
spread onto a Ni-foil. Two-electrode "coffee-bag" CugSns/Li-cells were made as described earlier Cell cycling was 
performed on a Digatron BTS-600 battery tester, in galvanostatic mode, with a current density of 0.022 mA/cm^. 

In situ X-ray diffraction data of CueSns electrodes in Li/CueSns cells were collected during the initial discharge 
in transmission mode using a STOE and CIE GmbH STADI powder diffractometer fitted with a position-sensitive 
detector (CuKai radiation). Measurements were recorded with a fixed detector covering 38.8° - 45.3° in 29 while a 
constant current of 0.15mA was applied by a MacPilell instrument. 

Table I: The optimized and the experimental |inl | atomic positions in the ?7'-Cu6Sn5 unit cell and the experimental 
atomic positions in the Li2CuSn unit cell The atomic designations are as in Ref. ^3 s-nd the site designations 

are in Wycoff notations J^] . 



Compound 


Atom 


site 


optimized positions 
X, y, z 


experimental positions 

X, y, z. 


ry'-CueSns 


Cul 


8f 


0.10200, 0.47500, 0.20610 


0.10096, 0.47297, 0.20236 




Cu2 


8f 


0.30570, 0.50600, 0.60650 


0.30620, 0.50404, 0.60972 




Cu3 


4a 


0, 0, 


0, 0, 




CuA 


4e 


0, 0.16800, i 


0, 0.16020, i 




Snl 


8f 


0.39000, 0.16500, 0.53000 


0.39106, 0.16250, 0.52867 




Sn2 


8f 


0.28000, 0.65000, 0.35650 


0.28518, 0.65499, 0.35792 




Sn3 


4e 


0, 0.80500, i 


0, 0.79892, i 


Li2CuSn 


Cu 


4c 




1 1 1 

4' 4' 4 




Sn 


4a 




0, 0, 




Li 


24f 




i 0, 




Li 


4c 




3 1 1 
4' 4' 4 



RESULTS AND DISCUSSION 

The structural optimization and the electric field gradients 

First, we shall discuss the importance of the atomic relaxations in r/'-CugSns. The base centered monoclinic 
structure of ?7'-Cu6Sn5 10] has four different kinds of Cu atoms at sites 4e, 4a and 8f and three different kinds of Sn 
atoms at site 4e and 8f with a total of 4 formula-units (44 atoms) per conventional unit cell (from now on referred 
to as the CUC). The calculations performed for the experimental atomic positions jl3 show that there are forces on 
the atoms which can not be neglected, so the relaxation of the structure is needed. On complete relaxation of the 
atomic positions, there is a gain in energy of 2.19 mRy/formula-unit of CugSns. The optimized atomic positions are 
presented in table I. 

The experimentally observed parameter that is sensitive to small structural changes and depends on the electronic 
charge distribution in a crystal is the electric field gradient (EFG). We have determined the EFG and the asymmetry 
parameter for both the relaxed and unrelaxed structures for comparison with future experimental work and the values 
are presented in table IL The details of the method of calculation of the EFG and the asymmetry parameter using 
FP-LAPW are given elsewhere [isl H^ . On relaxation of the structure the magnitude of Yzz for Snl and Sn2 
increases and decreases for Sn3. The asymmetry parameter shows an increase for Snl and Sn3 and decrease for Sn2. 
For all the four kinds of Cu atoms the magnitude of Yzz increases and the asymmetry parameter decreases on the 
structural relaxation of ry'-CugSns. 

In the case of Li2CuSn there are no forces on the atoms and so the atomic relaxations are not needed and all the 
further calculations are performed with the atoms at the ideal positions. Since the atoms in Li2CuSn are present at 
high symmetry positions of a cubic lattice the theoretical EFG vanishes. 

Table IL The electric field gradient (Vzz) and the asymmetry parameter (77) for different Cu and Sn sites in the 
ry'-CugSns unit cell with the optimized and the experimental |lOl] atomic positions. 
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Optimized structure 




Experimental structure 


[10] 


Atom 


Vzz 


(XlO^i V/m^) 


V 


Vzz (XlO^i V/m2) 




Cul 




-4.05 


0.326 


-3.84 


0.385 


Cu2 




-4.46 


0.153 


-4.32 


0.163 


Cu3 




-4.58 


0.257 


-4.40 


0.331 


CuA 




1.42 


0.508 


1.36 


0.707 


Snl 




-10.35 


0.743 


-10.10 


0.644 


Sn2 




4.18 


0.519 


3.80 


0.806 


Sn3 




-14.11 


0.718 


-15.02 


0.696 



The electron charge density 

Figure ?? and ?? show the difference between the crystalline charge and the superposed atomic charge for ?7'-Cu6Sn5 
in (102) plane and for Li2CuSn in (110) plane respectively. It is clear from figure 2 that in case of ?7'-Cu6Sn5 the Sn 
atoms gain a small fraction of the charge lost by the Cu atoms. It may be noted that for the Sn atoms the charge 
distribution in the inner atomic region (around the core) is spherical, while the increased charge in the outer atomic 
region is deformed. This deformation of the electron charge density is much more pronounced for the Cu atoms than 
the Sn atoms. The plot also shows the electron cloud between the atoms indicating metallic character of the bonds. 

On the other hand in the case of Li2CuSn, (figure 3 ) the Li atoms are strongly depicted of electrons while the 
Cu atoms show a gain in electronic charge. The nearly spherical shape of the charge distribution indicates an ionic 
character of the bonds. 

A comparison of the figures 2 and 3 indicates the modification in the bonding character. Lithiation of T^'-CueSns 
clearly results in a charge transfer from Li to the Cu and Sn atoms to form the more ionic Li2CuSn phase. In the 
battery, Li is inserted as an ion with the electrons supplied from the outer circuit redistributing in the structure. In 
this paper we have chosen to use the term Li atoms, since the calculations are based on Li atoms within the structures. 

The hthiation of r^'-CueSns 

According to the reactions (1) and (2), the total number of Li atoms reacting with one formula unit of CugSns is 
22. The experimental voltage vs. composition curve (figure 4) shows a maximum x of almost 22, which is in very 
good agreement with the suggested mechanisms. The in sz^m-XRD and Mossbauer measurements have, however, 
shown that the reactions taking place above 0.7V (section A in figure 4) can be attributed to surface reactions and 
reduction of oxides present in the sample due to the synthesis method (high-energy ball- milling) These reactions 
are irreversibly consuming Li and are not directly involved in reactions (1) and (2). The total experimental value for 
these reactions is thereby slightly smaller. The plateau, at around 0.4V (in section C of figure 4), is associated with 
the transformation of 77'-Cu6Sn5 to Li2CuSn. It has been suggested, that during this phase transformation, half of 
the Sn atoms in the ry'-CugSns structure are displaced into neighbouring Sn-strings 0,1^. The Li2CuSn-type phase 
at the end of the plateau is most likely lithium deficient, as further lithiation down to 0.2V (the sloping part in the 
voltage curve) results in a gradual expansion of the Li2CuSn cubic cell axis, i.e. a solid-solution behaviour within the 
structure. The process below 0.2V (section D in figure 4) involves lithiation of Sn to finally yield Li4.4Sn. We have 
focused mainly on the process described by reaction (1) (the CugSns- Li2CuSn transformation). 

From the total energy calculations the average insertion voltage (AIV) for the CufiSn.ts-Li^CuSn transformation is 
calculated. The details of the method of calculation of AIV are given elsewhere [la,lla|- The calculated value of the 
AIV is 0.378V which, is in good agreement with the experimental value of the plateau in the voltage profile (figure 4 
)• 

In the voltage profile, the region between 0.7V and 0.4V (section B in the figure 4), i.e. after the reduction of the 
oxides present in the sample and before the start of the plateau, there is a sloping region involving approximately two 
Li atoms. This region would then correspond to the first two Li atoms inserted in the Ty'-CugSns structure (i.e. x=l 
to x=2 in the reaction (1)). The in situ XRD measurements in this region show that the first lithium insertion does 
not result in any major structural rearrangements and the Ty'-CugSus peaks remain imaltered. The first-principles 
total energy and force calculations have been performed to study the various possibilities of how these two Li atoms 
are inserted into the Ty'-CugSns structure. In figure 5 the calculated total valence charge density of Ty'-CugSns in 
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the (102) plane is plotted (there exist four such planes per CUC). The minimum densities are found in the regions 
marked by A and B in the figure (there are 8 crystallographically symmetric sites of each kind per CUC). By looking 
at the charge density plots in various planes in the unit cell (not presented in the paper to keep the number of figures 
limited) the regions A, B and C are found to have minimum electron density and hence represents voids where the Li 
atoms can enter. 

With this information we go on to study how well these sites incorporates a Li atom. By total energy calculations 
it was found that the energy requirement for the insertion of a Li atom at site A in the CUC (one Li atom per CUC 
corresponds to x=0.25 in reaction (1) and x=3 in figure 4) is 41mRy/Li-atom. On the other hand the energy required 
for insertion of a Li atom at site B or C in the CUC is 66mRy/Li-atom and 120mRy/Li-atom respectively. This shows 
that it is energetically favourable for the first two Li atoms to enter at the sites A rather than the sites B or C. 

On complete relaxation of the structure with eight Li atoms per CUC (making a total of 52 atoms per CUC, i.e. 
x=2.0 in reaction (1) and x=4.5 or end of section B in figure 4) at the eight crystallographically symmetric A sites, two 
of which are marked in figure 5, the gain in energy is 190mRy per Li2Cu6Sn5 formula-unit with a volume expansion of 
4.6%. In this relaxation the Sn (Sn2) atoms move towards the site B, occupying the position with coordinates: 0.643, 
0.450, 0.200 and the Li atoms move to the position earlier occupied by the Sn (Sn2) atoms. This picture supports 
the earlier suggestions 0,0 about the structural rearrangement on lithiation of ry'-CugSns. However, with just two 
Li atoms per formula-unit there is no major structural rearrangement and the Sn atoms do not reach the sites B to 
form the Sn chains, but move towards sites B and come to rest nearly half way between B and its original position 
(Sn2). This rearrangement is not detected experimentally. 

There are two scenarios for further lithiation, one where lithium is inserted without any copper extrusion and the 
other where the lithium insertion is accompanied by extrusion of the CuA atoms (figure 5) from the structure. There 
exists no information about the morphology and the phase adopted by the extruded copper atoms, so the second case 
can not be studied theoretically. It was found that the first case (lithium insertion without copper extrusion i. e. 3 Li 
atoms per CugSns formula-unit) is accompanied by a volume expansion of 47.3 % and costs large amount of energy 
(290mRy/Li-atom) which in turn would imply a sloping voltage profile in the beginning of the section C i.e. between 
x=4.5 to x=5.5 in the figure 4 (which corresponds to x=2 to x=3 in reaction (1)). However, experimental voltage 
profile instead shows a plateau after insertion of 2 Li atoms per CugSns formula-unit. We, therefore, suggest that 
after the insertion of two Li atoms per CugSns formula unit, the Cu (CuA) atoms start to extrude from the structure 
accompanied by structural rearrangements. 

In situ XRD measurements were performed to closely follow the phase transformation from y/'-CueSns to Li2CuSn. 
The results can be seen in figure 6, showing diffraction patterns taken at consecutive steps at the 0.4V plateau {i.e. 
beyond 2 Li atoms per formula-unit). On lithiation, the CugSns peak at ~ 43° gradually decrease in intensity, while 
the Li2CuSn peak at ~ 40.5° starts to increase in intensity. The data shows that the plateau corresponds to the 
two-phase transformation between CugSns and Li2CuSn, where half of the Sn atoms move to the neighbouring Sn 
chains. This transformation is accompanied by quite large structural rearrangements. The CugSns peak decreases 
in intensity and broadens before the first sign of the Li2CuSn phase. The diffraction peak of the resulting Li2CuSn 
phase is even broader due to a more disordered phase. This supports the theoretical picture presented above. 

The Mossbauer and XRD H, ^3 experiments, performed for phase determination during lithiation of T^'-CugSns, 
are not sensitive to the point of copper extrusion and the phase adopted by the extruded copper atoms. Since the 
cycling stability depends upon the copper extrusion and stability of the structure we suggest EXAFS (Extended 
X-ray Absorption Fine-structure Spectroscopy) and HRTEM (High-Resolution Transmission Electron Microscopy) be 
performed to provide more information on the exact starting point of the copper extrusion as well as the structure 
and morphology of the extruded copper. Further, the knowledge of the phase of the extruded copper atoms would 
facilitate future theoretical work on the copper deficient lithiated yy'-CueSns. 

CONCLUSIONS 

Ty'-CugSns is a candidate anode material for lithium-ion batteries. At a discharge voltage above 0.2V lithium-ions 
are inserted in Ty'-CugSns reversibly to form Li2CuSn. In this region the cycling capacity is stable around 200 mAh/g. 

From the first principles total energy and force calculations the atomic positions in r^'-CugSns are optimised and 
EFG (the structure and electronic charge dependent parameter) is determined. The AIV of the Ty'-CugSns to Li2CuSn 
transformation has been calculated to be 0.378V which is in good agreement with the experimental result of ^ 0.4V. 

By calculating and plotting the difference between the valence- and atomic charge in a crystal plane, the change in 
the bonding nature on lithiation of ry'-CugSns could be studied. It was found that the ionicity of the bonds increases 
on going from r/'-CueSns (where the bonds have metallic character) to Li2CuSn. 
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Total energy and force calculations show that it is energetically favourable for the first two incoming Li atoms 
(two Li atoms per r^'-CueSns formula- unit) to enter at the crystallographically symmetric sites with the positional 
coordinate 0.907, 0.625, 0.0408. Complete relaxation of the structure leads to small structural rearrangement with 
the relaxation of the Sn (Sn2) atoms to occupy sites with the positional coordinate 0.643, 0.450, 0.200 and Li atoms 
to the positions earlier occupied by the Sn (Sn2) atoms. This leads to a large gain in energy of 190mRy/formula-unit 
of Li2Cu6Sn5 with a small volume expansion of 4.6%. Further addition of Li atoms (after two Li atoms per CueSns 
formula-unit) without any copper extrusion from the structure costs a large amount of energy leading to a sloping 
voltage profile. However, the experimental voltage profile shows a plateau after x=4.5 (which corresponds to x=2.0 
in reaction (1)), which suggests that after the insertion of two Li atoms per CueSns formula- unit (eight Li atoms per 
cue) the copper atoms arc extruded form the structure leading to large structural rearrangement. This picture is 
supported by our XRD measurements which does not show any detectable structural change up to 2 Li atoms per 
CueSns formula-unit but indicates a major structural rearrangement beyond that. 
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FIGURE CAPTIONS 

• Fig 1. The conventional unit cell of vy'-CugSns. The (102) plane is marked. 

• Fig 2. The difference between the crystalline valence charge and the superposed atomic charge for r^'-CueSns in 
the (102) plane. The difference is given in the units of charge per Oq, with the negative numbers in the legend 
indicating an increase in electrons. 

• Fig 3. The difference between the crystalline valence charge and the superposed atomic charge in the (110) 
plane of the Li2CuSn unit cell. The difference is given in the units of charge per Uq, with the negative numbers 
in the legend indicating an increase in electrons. 

• Fig 4. The voltage profile for the first discharge of a Li/CueSns cell, x indicates number of Li atoms inserted 
per CueSns formula unit. 

• Fig 5. The total valence electron charge density per Oq for Ty'-CueSns in the (102) plane. The six regions of low 
charge density are marked as A, B and C. 

• Fig 6. The in situ XRD data of a Li/CugSns cell measured at the 0.4V plateau during galvanostatic discharge. 
The Ni peak shown is from the current collector. 
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